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ABSTRACT 


Analytical formulations are presented which account for the coupled mechanical, electrical, and 
thermal response of piezoelectric composite laminates and plate structures. A layerwise theory is 
formulated with the inherent capability to explicitly model the active and sensory response of 
piezoelectric composite plates having arbitrary laminate configurations in thermal environments. 
Finite element equations are derived and implemented for a bilinear 4-noded plate element. 
Application cases demonstrate the capability to manage thermally induced bending and twisting 
deformations in symmetric and antisymmetric composite plates with piezoelectric actuators, and show 
the corresponding electrical response of distributed piezoelectric sensors. Finally, the resultant 
stresses in the thermal piezoelectric composite laminates are investigated. 


INTRODUCTION 


The development of piezoelectric composite materials offers great potential for use in advanced 
aerospace structural applications. By taking advantage of the direct and converse piezoelectric 
effects, piezoelectric composite structures can combine the traditional performance advantages of 
composite laminates along with the inherent capability of piezoelectric materials to adapt to their 
current environment. Extensive development of analytical methods for modeling the isothermal 
behavior of piezoelectric composite structures has been reported. The initial application of 
piezoelectric materials as actuators involved the vibration control of beams 1 ' 2 and lead to the 
development of simplified models for beams. Advances in theoretical models lead to both a 
piezoelectric plate theory 3 " 4 and a piezoelectric thin shell theory 5 , based on the assumptions of 
classical plate and shell theory, respectively. Subsequent developments led to a variety of finite 
element formulations for beam 6 " 7 , plate** 9 , shell 10 ' 13 , and solid 14 * 15 elements. 
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Additional approaches have been reported on refined piezoelectric laminate theories which account 
for the limitations of classical laminate plate theory (CLPT) 16 * 17 . These approaches, called discrete 
layer theories, provide a generalization of both the classical and higher order theories, and allow 
separate displacement fields to be assumed in each layer of the laminate. A discrete layer beam theory 
for analyzing piezoelectric beams was first reported by Robbins and Reddy 6 using an induced strain 
approach to approximate the piezoelectric effect. Heyliger and Saravanos 1 ** 20 subsequently presented 
discrete layer beam and plate theories explicitly accounting for the coupled equations of 
piezoelectricity, which have lead to generalized multi-field finite element formulations. 

All of the previously described works neglect the implication of thermal effects on both the active and 
sensory response of piezoelectric structures, even though an area where piezoelectric materials may 
provide dramatic advantages is in the development of smart thermal structures with the capability to 
sense and actively compensate for thermally induced deformations. Thermal loads will typically affect 
the response of smart piezoelectric laminates with four distinct physical mechanisms: (1) induction 
of thermal strains due to coefficient of thermal expansion mismatch between the various composite 
plies and piezoelectric layers, (2) pyroelectric effects on the electric displacement of the piezoelectric 
material, (3) changes in the piezoelectric and dielectric properties of the piezoelectric materials, and 
(4) changes in the elastic properties of the composite and piezoelectric materials. Only limited 
research has been reported in this area. T wo-dimensional thermopiezoelectric equations for plates 
were formulated by Mindlin 21 . A thermopiezoelectric laminate plate theory and a thin shell theory, 
based on classical assumptions, were reported by Tauchert 22 and Tzou and Howard 23 , respectively. 
Subsequent developments in thermopiezoelectric finite elements include a beam element by Rao and 
Sunar 3 ^ 25 , a shell element using induced strain approximations for the electric and thermal effects by 
Chandrashekhara and Kolli 26 , a solid element using an induced thermal strain approach by Ha et al. 27 , 
and a thin solid element by Tzou and Ye 2 ** 29 . Many of these previous approaches utilize simplified 
laminate assumptions and neglect the coupling effects existing in thermopiezoelectric laminates. To 
remedy these limitations, discrete layer approaches for smart thermopiezoelectric beam structures 
were reported by Lee and Saravanos 30 * 32 . 

This paper presents the extension of the discrete layer mechanics 30 * 32 for the analysis of smart 
thermopiezoelectric plate structures and addresses the problem of active thermal distortion 
management with smart piezoelectric plates. The mechanics account for the coupled mechanical, 
electrical, and thermal response of piezoelectric laminates at the material level through the thermo- 
piezoelectric constitutive equations. The displacements, electric potential and temperature are 
assumed to be layerwise (piecewise continuous) fields through the laminate thickness. This layerwise 
generalizatio n provides the capability to capture the locally induced piezoelectric effects, which leads 
to increased accuracy in laminate stress predictions (especially for thick laminates, and/or laminates 
with strong thermal and elastic inhomogeneities through-the-thickness). A corresponding finite 
dement formulation is presented using the layerwise laminate theory and a 4-noded plate element is 
developed. Numerical studies demonstrate the capability to actively manage thermally induced 
twisting and bending deformations in piezoelectric symmetric and antisymmetric composite plates 
subjected to thermal gradients, as well as examining the corresponding voltage response of 
piezoelectric sensors and the resulting intralaminar and interlaminar stresses in active thermo- 
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piezoelectric plates. 


GOVERNING MATERIAL EQUATIONS 

This section outlines the governing equations for thermopiezoelectric materials. The mechanical 
response is represented by the stress-equilibrium equation, 

+ ft - P (1) 

where o^, p, U* are the stress, body force per unit volume, density, and displacement, respectively; 
with i, j = 1, 2, 3. Maxwell's equation for the conservation of electric displacements D ; is 

D ui = 0 (2) 


The constitutive equations for a thermopiezoelectric material 33 employing standard contracted nota- 
tion are 
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where S, represents the strain tensor; E m the electric field vector; 0 = AT = T - T c is the temperature 
difference from the current temperature T and the thermally stress free reference temperature T c ; and 
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where C a(s and s ap are the elastic stiffness and compliance tensors; d am and e am are the different forms 
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of the piezoelectric tensor due to inversion; €„* is the electric permittivity tensor; a, and X a are the 
different forms of the coefficient of thermal expansion due to inversion; p m is the pyroelectric 
constant; superscripts E, a, S, and T, represent constant voltage, constant stress, constant strain, and 
constant temperature conditions, respectively; a, (3 = 1, ...,6; and k, m = 1, 2, 3. 

The small deformation strain-displacement relations are 

s v = \ <“« + V (*) 


and the electric field vector is related to the electric potential <f> by 


(9) 


Through use of the divergence theorem and neglecting body forces, Eqs. (1) and (2) can be expressed 
in an equivalent variational form as 


J (p u t bu t + o tj bS v - D ( 6E t ) dV = j t ( bu t dT + J q 6<J> dT 


( 10 ) 


where tj are the surface tractions applied on the bounding surface T; q is the electrical charge applied 
on the surface T p of the piezoelectric material; and V represents the whole volume including both 
composite and piezoelectric materials. 


DISCRETE LAYER LAMINATE THEORY 

A discrete layer laminate theory for thermopiezoelectric composite plates is formulated by 
introducing the following piecewise continuous approximations for the state variables 
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where N is the number of interpolation functions v|/ j . Currently, i|r j is represented using linear Lagran- 
gian interpolation functions, although in general any order of interpolation functions can be used. U\ 
V\ W°, (f^, and & are the generalized laminate state variables, in which a constant through-the- 
thickness displacement is assumed for w. This assumption of constant through-the-thickness 
displacement leads to a simplified set of constitutive equations, derived from Eqs. (3)-(6), shown in 
Appendix A 

By substituting Eqs. (5)-(9) and (1 1)-(15) into Eq. (10) and integrating through-the-thickness, the 
following variational form is obtained 
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in which the dependence of the z-coordinate has been separated into the generalized density matrix 
[P], the stiffness matrices [A], [B], [C], [D], the piezoelectric matrix [E], the dielectric permittivity 
matrix [G], and the thermal force matrices [f 5th ] and [q >th ]. These generalized laminate matrices are 
listed in Appendix B. 


FINITE ELEMENT FORMULATION 

The finite element formulation for a composite piezoelectric plate is obtained by incorporating 
additional local in-plane approximations to the state variables 

u(x,y,z,t) = £ £ U*{t) *'(x,y) +V) ( 17 ) 

<-i /-i 


6 


( 18 ) 


v(x ,y,z,f) = E E F*(f)*W) ♦>(*) 

/-I /-I 

M 

w(x,y,z,t) = E * - (0*W) 

/-I 


$(x,y,z,t) = E E ♦*(*) R '(^F) V(») 


<-l /-I 


(19) 

( 20 ) 


where M is the number of in-plane shape functions R. Currently, R is represented using bilinear 
Lagrangian interpolation functions. 


By combining Eqs. (16)-(20), the following finite element matrix formulation is obtained for the case 
of a constant through-the-thickness displacement plate 
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where the elements of the submatrices above are calculated in accordance with the generalized 
discrete layer laminate matrices defined previously in Appendix B. The finite element submatrices 
for the mass, stiffness, external force, and thermal force are included in Appendix C. 


The coupled finite element formulation can also be expressed in a compact form with the electric 
potential partitioned into active and sensory components such that 
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where the superscripts s and a indicate the partitioned electric potential vectors in either sensory or 
active configurations, respectively. This form has the advantage of positioning the unknown variables 
(displacements and sensory electric potentials) in the left-hand terms, while the known quantities 
(mechanical loads, thermal loads, electric charges, and active voltages) are included in the right-hand 
terms. 

The partitioned finite element formulation of Eq. (22) can be easily uncoupled into the following 
independent equations for the structural displacements 
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APPLICATIONS 


Simply Supported (0/±451 t Graphite/Epoxv Plate with Discrete Piezoceramic Patches 

This case study is a problem of active thermal distortion management, examined previously by Ha 
et al. 27 , of a 37.2 cm x 22.8 cm x 0.75 mm [0/±45] t graphite/epoxy plate with fifteen 6.0 cm x 6.0 
cm x 0.13 mm piezoceramic patches uniformly attached to each top and bottom surface (see 
Figure 1). A thermal gradient (50°C on the top surface and -50°C on the bottom surface) is 
applied to the plate. Figure 1 illustrates the geometry and finite dement mesh of the plate. The 
plate is simply supported along the two edges parallel to the y-axis and is free on the two edges 
parallel to the x-axis. The material properties used for both the composite and piezoceramic 
materials are listed in Table I. The objective was to apply increasing active voltages on both the 
top and bottom piezoceramic patches to minimize the out-of-plane deflection induced by the 
thermal gradient. Additional studies investigate the sensory response of the piezoceramic patches 
operating in a combined sensor and actuator mode, as well as the development of local stresses in 
the plate with all piezoceramic patches operating as actuators. 

Active Thermal Distortion Control . The resulting bending deformation with 0 Volts applied on 
both surfaces of the piezoceramic patches is shown in Figure 2(a). Through application of 
increasing active voltages in the piezoceramic patches, the thermally induced deformation of the 
plate ran be gradually eliminated. Figures 2(b) and 2(c) demonstrate the reduction in the out-of- 
plane deflection achieved by applying active voltages of 34 V and 70 V, respectively. The active 
voltage is applied on the outer surface of each of the piezoceramic patches, while the inner surface 
(i.e. the surface in contact with the plate) is grounded. The same voltages are applied to both the 
upper and lower piezoceramic patches. A comparison of centerline deflections (along y/b=0.5) 
predicted by Ha et al. 27 and by the current formulation, along with the corresponding active 
voltages, is shown in Figure 3. In order to obtain the same reduction in centerline deflection as 
Ha et al., slightly higher active voltages were applied using the current method (34 and 70 Volts, 
as compared to 31 and 61 Volts for Ha et al.). These differences between the two analyses are 
most lik ely due to the different formulations and elements used. Ha et al. uses a three-dimen- 
sional brick element in their analysis, while the current method utilizes a layerwise plate theory and 
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dement which assumes a constant through-the-thickness displacement w. This use of a plate theory 
neglects ce rtain material properties that Ha et al's. analysis considers (specifically Ci 3 , C^, C 33 , d u , 
d 33 , and a 33 ). Another difference is the explicit incorporation of thermal terms in the constitutive 
equations of the current formulation, which are neglected in Ha et al's. analysis. 

goncnry Vnitapfts Besides the fully active configurations previously studied (active potentials applied 
to all piezoceramic patches), other electric configurations yielding combinations of piezoceranuc 
sensors and actuators are sometimes more desirable. The current formulation has the capability to 
explicitly predict sensory voltages. These sensory capabilities are demonstrated on an active/ sensory 
configuration of the [0/±45], graphite/epoxy plate. The upper layer of piezoceramic patches are used 
as sensors, while the bottom layer of patches are maintained as actuators. In the sensory application, 
the inner surface of each piezoceramic patch is grounded while the outer surface is free to develop 
a sensory electric potential. Figures 4(a)-4(c) show the sensory voltages that develop on the upper 
piezoceramic patches upon combined application of a thermal gradient with active voltages of 0, 34, 
and 70 Volts on the lower piezoelectric patches, respectively. Generally, there is a reduction in the 
sensory voltage values which corresponds to the reduction in out-of-plane deflection with increasing 
active voltages. However, additional contributions remain from local stresses and pyroelectric 
effects. In actual applications, these measured sensory voltages will be used in a control loop to 
provide the inferred feedback essential for monitoring and actively controlling the thermal distortions. 

Local Stresses . The development of high local stress fields in a thermopiezoelectric laminate from 
the mismatch in induced thermal and piezoelectric strains is also an area of primary concern, since it 
will affect the integrity of these structures. As stated previously, one strong advantage of the 
layerwise approach is the accurate calculation of intralaminar in-plane stresses and interlaminar shear 
stresses in piezoelectric composite laminates (Heyliger and Saravanos 2 * > ). Consequently, this section 
presents the stress fields for the [0/±45], graphite/epoxy plate problem with both piezoceramic layers 
used as actuators. In these studies, a continuous piezoceramic layer is used instead of many 
distributed piezoceramic patches to eliminate the need for a highly refined mesh. The stresses are 
no ndimensi onalized using an equivalent laminate moduli defined as 

» 

<E Li}’ G Uj > = 7 / dz 

" 0 

Figures 5(a)-5(c) show the variation of normal stresses (o„, o w ) and in-plane shear stress (o^) 
through-the-thickness at x/a=0.5 and y/b=0.4 with the different applied voltages, while Figures 5(d)- 
5(e) show the variation of interlaminar shear stresses (o^ OyJ through-the-thickness near the free 
edge at x/a=0.9 and y/b=0.4. The normal stresses in Figures 5(a) and 5(b) increase as higher active 
voltages are applied reflecting the increase in piezoelectric strains, which indicates that one potential 
tradeoff for minimizing thermal distortions would be increased normal stresses. Thus, the stress state 
in thermopiezoelectric laminates (and particularly in the actuators) should be among the primary 
design considerations. The normal stresses can also be separated into individual components as 
shown in Figures 6(a) and 6(b) for the 70 V case (the shear stresses contain only an elastic component 
in the current study). The three components correspond to the three terms found on the right hand 
side of Eq. (5) and represent the individual stresses which are induced from the elastic, piezoelectric, 
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and thermal effects. This shows how the individual stress components are affecting the overall stress 
state of the plate. 


Cantilevered f45j/-45 2 l Graphite/Epoxv Plate with Discrete Piezoceramic Patches 

A numerical study was also conducted on a [45/453] antisymmetric graphite/epoxy laminate with 
attached piezoceramic patches. The plate has the same geometry and material properties used in the 
previous section and is clamped along the x=0 edge, while the other three edges remain free. A 
the rmal gradient (15°C on the top surface and - 15°C on the bottom surface) is applied to the plate 
with 0 Volts applied on both surfaces of the piezoceramic patches. The objective of this study is to 
minimize, both individually and in combination, the thermally induced bending and twisting 
deformations through application of active voltages. 

Active Thermal Distortion Control . Application of a 30°C thermal gradient, with all piezoceramic 
patches grounded, results in the combined bending and twisting deformation shown in Figure 7(a). 
Thro ugh application of active voltages of varying polarities in the piezoceramic patches, the bending 
and twisting deformations can be inhibited either individually or in combination. As demonstrated 
in the previous problem, the bending deformation can be managed through application of voltages 
which have the opposite polarity. By applying 40 V to both the top and bottom layers of piezo- 
ceramic patches, the bending deformation can be eliminated as shown in Figure 7(b). The twisting 
deformation remains unaffected and is now more apparent due to the removal of thermal bending. 
In a similar manner, the twisting deformation of the plate can be eliminated, without affecting the 
bending behavior as shown in Figure 7(c). This requires the application of voltages with the same 
polarities, specifically 145 V is applied to the top layer of piezoceramic patches, while -145 V is 
applied to the bottom layer of piezoceramic patches. Finally, in order to inhibit both the bending and 
twisting deformations simultaneously, the previous voltages are superposed. Thus, by applying 185 
V on the upper and -105 V on the lower piezoceramic patches, all deformations in the plate can be 
compensated, as shown in Figure 7(d). These results demonstrate the potential of piezoelectric 
structures to selectively manage desired thermal distortions. 

Sensory Voltages . The sensory capabilities of the current formulation are demonstrated on a 
combined active/sensory configuration of the [45/453] graphite/epoxy plate. The upper layer of 
piezoceramic patches are used as sensors, while the bottom layer of patches are maintained as 
actuators. Figure 8 shows the voltages that develop on the upper surface sensors upon combined 
application of a thermal gradient and 0 Volts on the lower piezoelectric actuators. These sensory 
voltages will provide the electric feedback required to monitor the state of deformation in the plate. 

Local Stresses . This section presents the stress fields for the active [45 3 /-45 3 ] graphite/epoxy plate 
problem. As before, a continuous piezoceramic layer is used instead of many distributed piezo- 
ceramic patches to eliminate the need for a highly refined mesh in this study. Figures 9(a)-9(c) show 
the variation of normal stresses (o B) o^) and in-plane shear stress (o v ) through-the-thickness at 
x/a=0.5 and y/b=0.4 with the various applied voltages, while Figures 9(d)-9(e) show the variation of 
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inteHamin ar shear stresses (o„, o^) through-the-thickness near the free edge at x/a=0.9 and y/b=0.4. 
The normal stresses in Figures 9(a) and 9(b) indicate that similar high stress states arise from 
minimizin g the twisting of the plate individually (145 V on top and -145 V on bottom layer) or 
simultaneously with bending (185 V on top and -105 V on bottom layer). In contrast, minimizing the 
bending behavior independently (40 V on both top and bottom layers) has little effect on the normal 
stresses. It is recalled that the twisting was compensated by taking advantage of the extension- 
twisting coupling in this laminate. The in-plane shear stress shown in Figure 9(c) indicates that 
regardless of the applied active voltage, high shear stresses will occur in the plate. The normal 
stresses can also be separated into the individual components as shown in Figures 10(a) and 10(b) 
for the case of combined active elimination of thermal twisting and bending (185 V applied on the 
upper piezoceramic patches and -105 V applied on the lower piezoceramic patches). The three 
components correspond to the three terms found on the right hand side of Eq. (5) and show the effect 
of the individual stresses on the overall stress state of the plate. 


SUMMARY 


Layerwise laminate and structural mechanics were described to model the coupled mechanical, 
electrical, and thermal behavior of smart piezoelectric composite laminates and plate structures. 
A finite element formulation and a bilinear 4-noded plate element were developed and encoded 
into prototype software. Thus, a necessary analytical and computational capability was developed 
to analyze the global and local response of smart piezoceramic plate structures operating in 
thermal environments. 

The accuracy of the formulation was compared with previously published analytical results. 
Numerical studies were performed on symmetric [0/±45] s and antisymmetric [45 3 /-45 3 ] composite 
plates subject to thermal gradient with discrete piezoceramic patches attached on both surfaces, 
which demonstrated the capability to actively manage thermally induced bending and twisting 
deformations. The corresponding electric voltages developed at distributed piezoelectric sensors 
were also calculated. Finally, the severe stress fields through-the-thickness of the composite 
plates were predicted and the contributions of the thermal and piezoelectric components 
quantified. The numerical studies have indicated the significance of thermal effects on the 
performance of piezoelectric structures in extreme temperature environments and demonstrated 
the capabilities of the mechanics to accurately model such behavior. 


11 



REFERENCES 


1 Bailey , T. and Hubbard, J.E., "Distributed Piezoelectric-Polymer Active Vibration Control of a 
Cantilever Beam," J. Guidance , Vol. 8, No. 5, 605-61 1, 1985. 

2 Crawley, E.F. and de Luis, J., "Use of Piezoelectric Actuators as Elements of Intelligent 
Structures," AIAA J., Vol. 25, No. 10, 1373-1385, 1987. 

3 Lee, C.K., "Theory of Laminated Piezoelectric Plates for the Design of Distributed Sen- 
sors/Actuators. Part I: Governing Equations and Reciprocal Relationships," J. Acoust. Soc. 
Amer., Vol 87, No. 3, 1144-1158, 1990. 

4 Wang, B.-T. and Rogers, C.A., "Laminate Plate Theory for Spatially Distributed Induced Strain 
Actuators," J. Comp. Mat., Vol. 25, 433-452, 1991. 

5 Tzou, H.S. and Gadre, M., "Theoretical Analysis of a Multi-Layered Thin Shell Coupled with 
Piezoelectric Shell Actuators for Distributed Vibration Controls," J. Sound Vib., Vol. 132, No. 3, 
433-450, 1989. 

6 Robbins, D.H. and Reddy, J.N., "Analysis of Piezoelectrically Actuated Beams Using a Layer- 
Wise Displacement Theory," Comput. Struct., Vol. 41, No. 2, 265-279, 1991. 

7 Shieh, RC., "Governing Equations and Finite Element Methods for Multiaxial Piezoelectric 
Beam Sensors/ Actuators," AIAA J., Vol. 32, No.6, 1250-1258, 1994. 

8 Hwang, W.-S. and Park, H.C., "Finite Element Modeling of Piezoelectric Sensors and 
Actuators," AIAA J., Vol. 31, No. 5, 930-937, 1993. 

9 Chandrashekhara, K. and Agarwal, AN., "Active Vibration Control of Laminated Composite 
Plates Using Piezoelectric Devices: A Finite Element Approach," J. Intel. Mat. Syst. Struct., Vol. 
4, 496-508, 1993. 

10 Lammering, R, "The Application of a Finite Shell Element for Composites Containing Piezo- 
electric Polymers in Vibration Control," Comput. Struct., Vol. 41, No.5, 1 101-1 109, 1991 

11 Tzou, H.S. and Tseng, C.I., "Distributed Dynamic Identification and Controls of Flexible 
Shells: Theory and Finite Element Development," Proceedings of the 3 1st Structures, Structural 
Dynamics, and Materials Conference, Long Beach, CA, 2265-2273, 1990. 

12 Tzou, H.S. and Tseng, C.I., "Distributed Piezoelectric Sensor/Actuator Design for Dynamic 
Measurement/Control of Distributed Parameter Systems: A Piezoelectric Finite Element 
Approach," J. Sound Vib., Vol. 138, No. 1, 17-34, 1990. 


12 



13 Tzou, H.S. and Ye, R, "Analysis of Laminated Piezoelectric Shell Systems with C° Piezoelectric 
Triangle Finite Elements," Proceedings of Adaptive Structures and Composite Materials: Analysis 
and Application, ASME, AD-Vol. 45, 113-124, 1994. 

14 Allik, H. and Hughes, T.J.R, "Finite Element Method for Piezoelectric Vibration," Int. J. Numer. 
Method Eng., Vol. 2, 151-157, 1970. 

15 Ha, S.K. and Chang, F.-K., "Finite Element Modeling of the Response of Laminated Composites 
with Distributed Piezoelectric Actuators," Proceedings of the 31st Structures, Structural Dynamics, 
and Materials Conference, Long Beach, CA, 2323-2330, 1990. 

16 Reddy, J.N., "A Generalization of Two-Dimensional Theories of Laminated Composite Plates," 
Commun. Appl. Numer. Methods, Vol. 3, 173-180, 1987. 

17 Reddy, J.N., "An Evaluation of Equivalent Single-Layer and Layerwise Theories of Composite 
Laminates," Compos. Struct., Vol. 25, 21-35, 1993. 

18 Heyliger, P.R and Saravanos, D.A., "On Discrete-Layer Mechanics for Health-Monitoring 
Applications in Sensory/ Active Composite Laminates," Proceedings of ASME Winter Annual 
Meeting, New Orleans, LA, 1993. 

19 Heyliger, P.R., Ramirez, G., and Saravanos, D.A., "Coupled Discrete-Layer Finite Elements for 
Laminated Piezoelectric Plates," Commun. Numer. Methods Eng., Vol. 10, 971-981, 1994. 

20 Saravanos, D. A. and Heyliger, P.R, "Layerwise Mechanics and Finite Element for the Dynamic 
Analysis of Piezoelectric Composite Plates," Int. J. Solids Struct., in press. 

21 Mindlin, RD., "Equations of High Frequency Vibrations of Thermopiezoelectric Crystal Plates," 
Int. J. Solids Struct., Vol. 10, 625-632, 1974. 

22 Tauchert, T.R., "Piezothermoelastic Behavior of a Laminated Plate," J. Thermal Stresses, Vol. 
15, 25-37, 1992. 

23 Tzou, H.S. and Howard, R.V., "A Piezothermoelastic Thin Shell Theory Applied to Active 
Structures," J. Vib. Acoust., Vol. 116, 295-302 , 1994. 

24 Sunar, M. and Rao, S.S., "Distributed Thermopiezoelectric Sensors and Actuators in Structural 
Design," Proceedings of the 33rd Structures, Structural Dynamics, and Materials Conference, 
Dallas, TX, 890-895, 1992. 

25 Rao, S.S. and Sunar, M., "Analysis of Distributed Thermopiezoelectric Sensors and Actuators in 
Advanced Intelligent Structures," AIAA J., Vol. 31, No. 7, 1280-1286, 1993. 


13 



26 Chandrashekhara, K. and Kolli, M., "Thermally Induced Vibration of Adaptive Doubly Curved 
Composite Shells with Piezoelectric Devices," Proceedings of the 36th Structures, Structural 
Dynamics, and Materials Conference , New Orleans, LA, 1628-1636, 1995. 

27 Ha, S.K., Keilers, C., and Chang, F.-K., "Finite Element Analysis of Composite Structures 
Containing Distributed Piezoceramic Sensors and Actuators," AIAA J., Vol. 30, No. 3, 772-780, 
1992. 

28 Tzou, H.S. and Ye, R., "Piezothermoelasticity and Precision Control of Piezoelectric Systems: 
Theory and Finite Element Analysis," J. Vib. Acoust., Vol. 116, 489-495, 1994. 

29 Tzou, H.S. and Ye, R., "Pyroelectric and Thermal Strain Effects of Piezoelectric (PVDF and PZT) 
Devices," Proceedings of Adaptive Structures and Composite Materials: Analysis and Application, 
ASME, AD- Vol. 45, 125-132, 1994. 

30 Lee, H.-J. and Saravanos, D. A., "On The Response of Smart Piezoelectric Composite Structures 
in Thermal Environments," Proceedings of the 36th Structures, Structural Dynamics, and Materials 
Conference, New Orleans, LA, pp. 2876-2885, 1995. 

31 Lee, H.-J. and Saravanos, D.A., "A Coupled Layerwise Analysis of the Thermopiezoelectric 
Response of Smart Composite Beams," NASA Technical Memorandum 106889, 1995. 

32 Lee, H.-J. and Saravanos, D.A., "Coupled Layerwise Analysis of Thermopiezoelectric Composite 
Beams," AIAA J. , in press. 

33 Nye, J.F., Physical Properties of Crystals, The Clarendon Press, Oxford, 1964. 


14 



TABLE I: Material properties of piezoceramic (PZT G1195N) and graphite/epoxy 

(T300/976) composite 27 



Piezoceramic 

Graphite/Epoxy 

Elastic Moduli (GPa): 

E u 

63.0 

150.0 


63.0 

9.0 

E33 

63.0 

9.0 

Poisson's Ratio: 

U 12 

0.3 

0.3 

U 23 

0.3 

0.3 

«31 

0.3 

0.3 

Shear Moduli (GPa): 

G 12 

24.2 

7.10 

Gj3 

24.2 

2.50 

G31 

24.2 

7.10 

Density (kg/m 3 ): 

P 

7600. 

1600. 

Piezoelectric Charge Constant (pm/V) : 

d» 

0. 

0. 

^15 

0. 

0. 

d 24 

0. 

0. 

^25 

0. 

0. 

^31 

254. 

0. 

^32 

254. 

0. 

^36 

0. 

0. 

Electric Permittivity (nf/m): 

e n 

15.3 

0. 

e 22 

15.3 

0. 

€ 33 

15.0 

0. 

Thermal Expansion Coefficient (pm/m 

«n 

°C): 

0.9 

1.1 

a 22 

0.9 

25.2 

Reference Temperature, T„, (°C): 

20. 

20. 
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Appendix A: 


Thermopiezoelectric Constitutive Equations for a Constant Through-the-Thickness 

Displacement Plate 

The simplified constitutive equations for a constant through-the-thickness displacement plate are 




(A-2) 




(A-4) 
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Appendix B: 

Gene raliz ed Laminate Matrices for a Constant Through-the-Thickness Displacement Plate 


If L is the number of composite plies and piezoelectric layers, the laminate density matrix is 

P* = £ / P ♦*(*) i T(z) dz (B-l) 


/-x 


for ij = 1 1 and 22, while 


^33 = £ f P dZ 

l-l j 


The laminate stifl&iess matrices are 

A? - £ / C„ ♦‘(z) f(r) * 

/-I 


for ij = 1 1, 12, 16, 22, 26, 44, and 66; 

B k = y f c a * t(z) 

* ~ J « dz 


(B-2) 


(B-3) 


(B-4) 


- £ / c, * 

J-l , 


O*" = f [ C 5]T(£) ^ 

^ J • dz dz 


l-l 


for ij = 44, 45, and 55. The laminate piezoelectric matrices are 

c - i / •, ♦*<*> * 


(B-S) 


(B-6) 


(B-7) 


, km 


i-i x 


♦-(z) * 
5r 


(B-8) 


for ij = 3 1, 32, and 36. The laminate dielectric permittivity matrix is 
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(B-9) 



/ +‘(z) ♦"’(z) dx 


for ij = 1 1 and 22, while 



d»*0) d»"(z) ^ 
c>z dz 


The laminate thermal force vector is 

>& “ E / K ♦/ e , ♦*(*> 

>i i 

for i = 1 and 2. The laminate thermal electric displacement vector is 

- £ / />, */ e > +*(*) * 

w x 


for i =1 and 2, while 



(B-10) 


(B-ll) 


(B-12) 


(B-13) 
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Appendix C: 

Finite Element Submatrices for a Constant Through-the-Thickness Displacement Plate 


The mass submatrices are 


[M tj T = / [P v ] R F (x,y) R 9 (x,y) dA 


The stiffness submatrices are 


[* n r ■ / ( m„i 


dR p dR 9 


* [A 16 \ ( 


3R P dR 9 dR p dR 9 


dx dx 
* [Z>„] R'R* * MJ 

dR p dR 1 


dx dy dy dx 
**' ><M 


IV - / < M J ^ ^ ^ * UJ 


dx dy 
* W 45 ] R p R 9 + [AJ 


dy 

dy 

dR * 

dR 9 

dx 

dx 

dR p 

dR 9 

dy 

dx 


dR” dR l 
dy dy 


} dA 


iK»r = / { [^ 4S i RP + lB ” ] R 
* & 


a*: * r»_.i *1 

dx 


} aa 


iK„r - / ( [£„] ^- R< + 1 £ »] ^- R,)dA 


W ft MJ ^ ^ ^ < 


dy dy 


as* aR^ Mi Mi 

dx dy dy dx 


) 


* to,,]*’*' * l A «> *ZT > 

dc dx 


(C-l) 

(C-2) 

(C-3) 

(C-4) 

(C-5) 

(C-6) 
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(C-7) 


i*»r = / ( [^i *' + 

i ^ ar 


i* 24 r = / { [£ 32 ] ** + ^ * f > ^ 

A dy dx 


— dR p dR 9 r — ,a*' a** a*' a« ? 


IV - / u!& • KJ < 


dy dy dx dy dy dx 

dR ' dR 9 


) 


+ [C.J 


0* 0X 


> dA 


dR r dR q 


dR p dR q 


[*„]” ■ / < [<?„] *=- "T” * [OJ ^ *T“ * IO M ] R’R')dA 


dx dx 


dy dy 


The external forces are 


[FJ = / t t R p (x,y) dA 


The thermal forces are 


[QJ - / ( 9uk 


* / q R p (*,y) dA 

A , 

II 

a*'^) 

dx 

- //» 
4 

dR p {x iV ) 

dy 

0U' 

dR p 

ex * * 


(C-8) 

(C-9) 

(C-10) 

(C-ll) 

(C-12) 

(C-13) 

(C-14) 

(C-15) 
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Figure 1: Model of Graphite/Epoxy Plate with Attached Piezoceramic Patches 

(a) Top View 



Figure 1: Model of Graphite/Epoxy Plate with Attached Piezoceramic Patches 

(b) Front View 
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Figure 2: Active Compensation of Thermal Deformation of [0/±45] s Plate Subject to 
100°C Gradient (a): 0 Volts Applied on Upper and Lower Piezoceramic Patches 
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Figure 2: Compensation of Thermal Deformation of [0/±45] s Plate Subject to 100°C 
Gradient (b): 34 Volts Applied on Upper and Lower Piezoceramic Patches 
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Figure 2: Compensation of Thermal Deformation of [0/±45] s Plate Subject to 100°C 
Gradient (c): 70 Volts Applied on Upper and Lower Piezoceramic Patches 


24 



max 


Ha et al.: • V 

• Currents • V 

Ha et al.: 31 V 

A Current: 34 V 
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Ml 0.20 0.40 0.60 0.80 1.00 

x/a 

Figure 3: Centerline (y/b=0.5) Deflections of [0/±45] s Plate 
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Figure 4: Sensory Voltages on Upper Piezoceramic Patches of a [0/±45] s Plate 
Subject to 100°C Gradient (a): 0 Volts Applied on Lower Piezoceramic Patches 
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Figure 4: Sensory Voltages on Upper Piezoceramic Patches of a [0/±45] s Plate 
Subject to 100°C Gradient (b): 34 Volts Applied on Lower Piezoceramic Patches 
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Figure 4: Sensory Voltages on Upper Piezoceramic Patches of a [0/±45] $ Plate 
Subject to 100°C Gradient (c): 70 Volts Applied on Lower Piezoceramic Patches 
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z/h 

Figure 5: Through-The-Thickness Stresses in [0/±45] s Plate Subject to 100°C 
Gradient with Different Applied Voltages on Upper and Lower Piezoceramic 

Patches (a): 
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Figure 5: Through-The-Thickness Stresses in [0/±45] s Plate Subject to 100°C 
Gradient with Different Applied Voltages on Upper and Lower Piezoceramic 

Patches (b): 
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Figure 5: Through-The-Thickness Stresses in [0/±45] s Plate Subject to 100°C 
Gradient with Different Applied Voltages on Upper and Lower Piezoceramic 

Patches (c): a v 
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Figure 5: Through-The-Thickness Stresses in [0/±45] s Plate Subject to 100°C 
Gradient with Different Applied Voltages on Upper and Lower Piezoceramic 

Patches (d): o^ 
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Figure 5: Through-The-Thickness Stresses in [0/±45] s Plate Subject to 100° C 
Gradient with Different Applied Voltages on Upper and Lower Piezoceramic 

Patches (e): o^ 


33 





Total 

Elastic 

Thermal 

Piezoelectric 


z/h 


Figure 6: Elastic, Thermal, and Piezoelectric Stress Components in [0/±45] s Plate 
Subject to 100°C Gradient with 70 Volts Applied on Upper and Lower Surface 

Piezoceramic Patches (a): o ffl 
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Figure 6: Elastic, Thermal, and Piezoelectric Stress Components in [0/±45] s Plate 
Subject to 100° C Gradient with 70 Volts Applied on Upper and Lower Surface 

Piezoceramic Patches (b): o w 
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Figure 7: Active Compensation of Thermal Bending and Twisting Deformation of 
[ 45 3 /_ 45 3 ] plate Subject to 30 °C Gradient (a): Initial Thermal Deformation: 0 Volts 
Applied on Upper and Lower Piezoceramic Actuators 


36 



Figure 7: Active Compensation of Thermal Bending and Twisting Deformation of 
[ 45 3 /- 45 3 ] plate Subject to 30°C Gradient (b): Active Compensation of Thermal 
Bending: 40 Volts Applied on Upper and Lower Piezoceramic Actuators 
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Figure 7: Active Compensation of Thermal Bending and Twisting Deformation of 
[45 3 /-45 3 ] Plate Subject to 30°C Gradient (c): Active Compensation of Thermal 
Twisting: 145 Volts Applied on Upper and -145 Volts on Lower Piezoceramic 

Actuators 


38 



Figure 7: Active Compensation of Thermal Bending and Twisting Deformation of 
[45 3 /_45 3 ] Plate Subject to 30 °C Gradient (d): Complete Compensation of Thermal 
Deformation: 185 Volts Applied on Upper and -105 Volts on Lower Piezoceramic 

Actuators 
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Figure 8: Sensory Voltages on Upper Piezoceramic Patches of a [45 3 /-45 3 ] Plate 
with 0 Volts Applied on Lower Piezoceramic Patches 
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Figure 9: Through-the-Thickness Stress in [45 3 /-45 3 ] Plate Subject to 30°C Gradient 
with Different Applied Voltages on Upper and Lower Piezoceramic Patches (a): 
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Figure 9: Through-the-Thickness Stress in [45 3 /-45 3 ] Plate Subject to 30°C Gradient 
with Different Applied Voltages on Upper and Lower Piezoceramic Patches (b): 
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Figure 9: Through-the-Thickness Stress in [45 3 /-45 3 ] Plate Subject to 30 °C Gradient 
with Different Applied Voltages on Upper and Lower Piezoceramic Patches (c): 
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Figure 9: Through-the-Thickness Stress in [45 3 /-45 3 ] Plate Subject to 30 °C Gradient 
with Different Applied Voltages on Upper and Lower Piezoceramic Patches (d): 
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Figure 9: Through-the-Thickness Stress in [45 3 /-45 3 ] Plate Subject to 30 °C Gradient 
with Different Applied Voltages on Upper and Lower Piezoceramic Patches (e): o^ 
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Figure 10: Elastic, Thermal, and Piezoelectric Stress Components in [45 3 /-45 3 ] Plate 
Subject to 30°C Gradient with 185 Volts Applied on Upper and -105 Volts Applied 

on Lower Piezoceramic Actuators (a): 
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Figure 10: Elastic, Thermal, and Piezoelectric Stress Components in [45 3 /-45 3 ] Plate 
Subject to 30 °C Gradient with 185 Volts Applied on Upper and -105 Volts Applied 

on Lower Piezoceramic Actuators (b): o >7 
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